Spinning-particle model for the Dirac equation and the relativistic Zitterbewegung 
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We construct the relativistic particle model without Grassmann variables which meets the fol- 
lowing requirements. A) Canonical quantization of the model implies the Dirac equation. B) The 
variable which experiences Zitterbewegung, represents a gauge non-invariant variable in our model. 
Hence our particle does not experience the undesirable Zitterbewegung. C) In the non-relativistic 
limit spin is described by three-vector, as it could be expected. 
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I. INTRODUCTION 



The Dirac spinor ^ can be used to construct the 
four-dimensional current vector, ^r^^, which preserves 
for solutions to the Dirac equation, d^(^/r^^) = 0. 
As a consequence, the integral over space of its null- 
component, J cPx^^, does non depend on time. Hence 
the quantity ty'ty > admits the probabilistic interpreta- 
tion, and we expect that one-particle sector of the Dirac 
equation can be described in the framework of relativistic 
quantum mechanics. Although the true understanding of 
spin is achieved in the framework of quantum electrody- 
namics, a lot of efforts has been spent in attempts to 
construct semiclassical description of relativistic spin on 
the base of mechanical models [1-12]. However, the rel- 
ativistic quantum mechanics applied to the one-particle 
solutions predicts some controversial effects, like the Zit- 
terbewegung of the center-of-charge position operator [1, 
2] and the Klein's paradox. 

Analyzing the applicability of quantum-mechanical 
treatment to the free Dirac equation, Schrodinger no- 
ticed [1] that the center-of-charge position operator in the 
Heisenberg picture experiences rapid oscillations called 
Zitterbewegung. If we take the state vector with positive 
and negative energy components, the expectation value 
of the operator has a similar behavior. It is often assumed 
that Zitterbewegung represents the physically observable 
motion of a real particle [18, 19]. The analogous systems 
that are described by a Dirac-type equation and simu- 
late Zitterbewegung are under intensive study in different 
physical set-ups, including graphene, trapped ions, pho- 
tonic lattices and ultracold atoms (see [22] and the ref- 
erences therein). In this work we show that the status 
of this phenomenon in relativistic quantum mechanics is 
not as clear as is commonly believed. 

Besides the center of charge, x, in the Dirac theory we 
can construct the center-of-mass (Pryce-Newton-Wigner) 
[3, 4] operator x in such a way that the conjugated mo- 
mentum of x turns out to be the mechanical momentum 
for x. So the Dirac particle looks like a kind of composed 
system (this picture has been used by Schrodinger [1] to 



identify spin with inner angular momentum of the sys- 
tem). It further complicates the semiclassical analysis, as 
the Dirac equation gives no evidence as regards which of 
these two operators should be identified with the position 
of the particle. 

To understand the controversial properties of the one- 
particle Dirac equation, it would be desirable to have at 
our disposal the semiclassical particle model which leads 
to the Dirac equation in the course of canonical quanti- 
zation. It implies, in particular, that T -matrices should 
be produced through quantization of some set of classi- 
cal variables. While the problem has a long history (see 
[3-12] and the references therein), there appears to be no 
wholly satisfactory solution to date. The main difficulty 
consists of the proper choice of the basic classical vari- 
ables for construction of the spin space. If we start from 
some classical-mechaniccs action functional, the phase- 
space variables, say u a , 7Tb, necessarily obey the Poisson 
bracket algebra {oj a , Kb} = <5 a &- The number of variables 
and the algebra are different from the number of spin 
operators and their commutators (for instance, for non- 
relativistic spin they are [Si, Sj] = iheijkSk). To improve 
this, we need to impose constraints as well as to pass from 
the initial to some composed variables. This implies the 
use of the Dirac machinery for constrained theories [13- 
16]. Although the non-relativistic spin can be described 
along these lines [17], it seems to be surprisingly difficult 1 
to construct, in a systematic way, a consistent model for 
the relativistic spin. 

In this letter we propose the relativistic particle model 
without Grassmann (anticommutative) variables which 
leads to the Dirac equation. We apply our model to 
analysis of the relativistic Zitterbewegung, presenting a 
simple semiclassical argument on the non-physical char- 
acter of this phenomenon, which in our model represents 
the dynamics of an unobservable variable. 



'Electronic address: alexei.deriglazov@ufjf.edu.br; On leave of ab- 
sence from Dep. Math. Phys., Tomsk Polytechnical University, 
Tomsk, Russia. 



1 The Berezin-Marinov model [9] is based on the Grassmann vari- 
ables and leads to the Dirac equation. The problem here is that 
the Grassmann classical mechanics represents a rather formal 
mathematical construction. It leads to certain difficulties [9, 14] 
in attempts to use it for description of the spin effects on the 
semiclassical level, before the quantization. The Hanson-Regge 
theory [8] does not produce T-matrices. The Barut-Zanghi model 
[10] does not imply the Dirac equation. 
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II. MODEL-INDEPENDENT CONSTRUCTION 
OF THE RELATIVISTIC SPIN SURFACE 

To construct the relativistic spin surface, we start from 
the Dirac equation 



71 BCjAD + lf u J 



BD jAC\ 



r"(p^ + -A^+mc 



* = 0. 



Applying the operator T^(p^ + -A^) — mc, this implies 
the Klein- Gordon equation with non-minimal interaction 



e ,„. 2 eh 



(0" + -A»Y + ^F^T^ + m 2 c 2 
c 2c 



* = 0, 



where 



(2) 



(3) 



and Ffj, v = d^Ay—dyA^. We use the representation with 
hermitean T° and antihermitean T l 



1 

-1 



a 1 
-a 1 



(4) 



then [r^,r ,y ]+ = -2rf", rf v = (- + ++), and r°P, T° 
are the Dirac matrices a\ (3 [2]. We take the classical 
counterparts of the operators i M and p^ — —ihd^ in the 
standard way, which are p u , with the Poisson brackets 
{x»,p v } = -rf v . 

Let us discuss the classical variables that could pro- 
duce the r -matrices. The relativistic equation for the 
spin precession is usually obtained including the three- 
dimensional spin vector S % either into the Frenkel tensor 
<§>nv ; ^ u Uv = 0, or into the Bargmann-Michel-Telegdi 
four-vector 2 5 M , S^u^ = 0, where u v represents four- 
velocity of the particle. However, the semiclassical mod- 
els based on these schemes do not lead to a reasonable 
quantum theory, as they do not produce the Dirac equa- 
tion through canonical quantization. We now present 
arguments as to how this can be achieved in the formula- 
tion that implies inclusion of S % into the 5*0(2, 3) angu- 
lar momentum tensor J AB of five-dimensional space, A = 
(/Li, 5) = (0, 1,2,3,5), with the metric rj AB = (- + + + -). 

First, we analyze commutators of the T -matrices. The 
commutators do not form closed Lie algebra, but produce 
5*0(1, 3) -Lorentz generators ©. The set T M , forms 
closed algebra 

pr"", r^] = 2i(r 1 f " a T 1 ' 13 - v ^r ua - rf a v^ + T^ 5 r" a ).(5) 

The algebra can be identified with the five-dimensional 
Lorentz algebra 5*0(2,3) with generators J AB 

{ jAB^ j CD] = 2i{ n AC J BD - r, AD J BC - 



(6) 



assuming = J 5 ^, r"" = P v . 

To reach the algebra starting from a classical- 
mechanics model, we introduce ten-dimensional "phase" 
space of the spin degrees of freedom, uj a , n B , equipped 
(1) with the Poisson bracket {oj a , ir B } — i] AB . Consider the 



inner angular momentum 



J' 



EE 2{LO A TT B - U B TT A ) 



(7) 



Poisson brackets of these quantities form the algebra 

tACjBD _ ^AD jBC _ 
\BCjAD +rj BDjAC ) (g) 



{J AB ,J CD } 



PB 



Comparing © with ^ we conclude that the operators 
p^i ^ p/iw could be obtained by quantization of J AB . 

Since J AB are the variables which we are interested in, 
we try to take them as coordinates of the space uj a ,tt b . 
The Jacobian of the transformation (lu a , tt b ) — > J AB has 
rank equal seven 3 . So, only seven among ten functions 
J AB (to, 7r), A < B, are independent quantities. They can 
be separated as follows. By construction, the quantities 
obey the identity & va ^ J 5 U J al3 = 0, this can be solved 
as 



jij = (j50yl(j5ij0j _ jSjjOij 



(9) 



Hence we can take J 5m , J° l as the independent variables. 
We could complete the set up to a base of the phase space 
(lo a , tt b ) adding three more coordinates, for instance w 3 , 
a; 5 , 7T 5 . Quantizing the complete set we obtain, besides 
the desired operators J 5m , J°\ some extra operators w 3 , 
w 5 , 7T 5 . They are not present in the Dirac theory, and 
are not necessary for description of spin. So we need 
to reduce the dimension of our space from ten to seven 
imposing three constraints. There is one important re- 
striction on the choice of constraints. Canonical quanti- 
zation of a system with constraints implies replacement 
the Poisson by the Dirac bracket, the latter is constructed 
with help of the constraints. We need 50(2, 3) -invariant 
constraints T a , {T a , J ab }pb — 0, otherwise the Dirac- 
brackets algebra will not coincide with those of the Pois- 
son, ©. 

The only quadratic 50(2, 3) -invariants which can be 
constructed from uj a , n B are uj a u>a, ^> A t^a and tt a tta- 
So we restrict our model to live on the surface defined by 
the equations 



T 3 = tt a wa + a 3 



0; 



(10) 



T 4 = uj a lu a + a 4 = 0, T 5 =uj A ir A =Q, (11) 
where 03, 04 are some numbers. 



2 The conditions <&v- v u v = and S^u^ = guarantee that in 
the rest frame only three components of these quantities survive, 3 The rank has been computed using the program: Wolfram Math- 
which implies the right non-relativistic limit. ematica 8. 
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The matrix x! a'^bT 5 '" ^ has rank equal ten. So 



the quantities 



J 5 ", J°\ T 4 , T 5 , W 5 , 



(12) 



can be taken as coordinates of the space (uj a , ir B ). The 
equation J AB = 2(oj a tt b — uj b tt a ) implies the identity 

J ab Jab = 8[(^) 2 (tt b ) 2 - (u; A n A ) 2 } = 

8[(T 4 -a 4 )(T 3 -a 3 )~(T 5 ) 2 }, (13) 

then the constraint T3 can be written in the coordinates 
(H2ll as follows: 



n = 



(J AB ) 2 + 8(T 5 )' 
8(T 4 - a 4 ) 



03, 



(14) 



where J y are given by Eq. ([9]). Note that T3 does not 
depend on lu 5 . On the hyperplane T4 = T5 = it reduces 
to 



- 8a 4 T 3 = {J AB f - 8a 3 a 4 = 



(15) 



Eq. (p~5)) states that the value of SO(2, 3) -Casimir oper- 
ator 4 (J AB ) 2 is equal to 8a 3 a 4 . 

In the dynamical model constructed below, the equa- 
tion T3 — appears as the first-class constraint. It im- 
plies that we deal with a theory with local symmetry, 
with the constraint being the generator of the symmetry 
[20] . The coordinate w 5 is not inert under the symmetry, 
Suj 5 ~ {T3,oj 5 } ^ 0. Hence ui 5 is gauge non-invariant 
variable. 

Summing up, we have restricted dynamics of spin 
on the surface (fT0|) . (fTTj) . If (fT2|) are taken as coordi- 
nates of the phase space, the surface is the hyperplane 
T 4 = T 5 = with the coordinates J 5fJ -, J°\uj 5 subject to 
the condition (|15[) . Since oj 5 is gauge non-invariant coor- 
dinate, we can discard it. It implies that we can quantize 
J 5 ^, J° l instead of the initial variables lu a , tt b . 

Following the canonical quantization paradigm, the 
variables must be replaced by Hermitian operators 5 with 
commutators resembling the Poisson bracket 



[,]=iK{,}\ 



(16) 



Similarly to the case of T -matrices, brackets of the vari- 
ables J 5ai , J° 4 do not form closed Lie algebra. The non 
closed brackets are 



{J 5i , J 53 } = { J°\ J 01 } = -2JV, 



(17) 



where J lJ is given by Eq. © . Adding them to the initial 
variables, we obtain the set J AB — (J 5a \ J Ql , J u ) which 
obeys the desired algebra ©. 



According to Eqs. ([6]), <j8j) the quantization is achieved 
replacing the classical variables J 5/i , J^ v on T -matrices 6 . 
We assume that u A has a dimension of length, then J AB 
has the dimension of the Planck's constant. Hence the 
quantization rule is 



(18) 



This implies that the Dirac equation ([TJ can be pro- 
duced by the constraint (since the Zitterbewegung is a 
property of the free Dirac equation, we take temporarily 



A. 



0) 



To 



J 5/i + mch = 0. 



(19) 



In quantum theory, the Dirac equation implies the Klein- 
Gordon one. In contrast, in the classical theory the con- 
straint (TlTH) does not imply the mass-shell condition 



Ti = p 2 + m 2 c 2 = 0. 



(20) 



To improve this, we are forced to look for a classical 
model that produces this equation as an independent 
constraint. The model without the constraint (|20l) has 
been considered in [21]. It shows the same undesirable 
properties as those of Dirac equation in the classical limit. 



III. SPINNING-PARTICLE ACTION, 
CANONICAL QUANTIZATION AND THE 
DIRAC EQUATION 

According to the previous section, to describe the rel- 
ativistic spin we need a theory that implies the Dirac 
constraints (fTOj) . (fTTj) . (IT9 | . ([20)) . 

We recall that the Hamiltonian action for a sys- 
tem with the phase-space variables Q a , P a reads S = 
P a Q a — H, H = Hq + X a ^a, where Hq is the canoni- 
cal Hamiltonian and A a are the Lagrangian multipliers 
for the primary constraints $ Q . For the present case, we 
propose to consider the total Hamiltonian of the form 
H = \e a T a + A 

CO 

7r ea , where 7r ea are conjugate momenta 
for the auxiliary variables e a , the latter are associated 
with the constraints T a , a = 1, 2, 3, 4. So, let us consider 
the following d — 4 Poincare invariant action 



S 



drp^x 



^(p 2 +m 2 c 2 ) 



mch) + yT 3 + yT 4 + 7T eQ (A e 



e«) ,(21) 



If we omit the spin-space coordinates, ui A = n A = 0, 
Eq. ([2"Tj) reduces to the well known action of a spinless 
rclativistic particle 

So = I drp^ 1 - \e{p 2 + m 2 c 2 ) + n e (e - A e ). (22) 



4 In quantum theory, for the operators (0 we have: 
J AB Jab = 20h 2 . 

5 The matrices , are Hermitian operators with respect to 6 Replacing ([9} by an operator J IJ ' (T^ , F° l ) we arrange the oper- 
the scalar product *2) = Vpjr ^. ators F in such a way, that J ij (r) = r ij . 
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Variation of the action ([21]) with respect to e a leads 
to the constraints T a = 0. Preservation in time of the 
constraint T 4 , T 4 = 0, implies T5 = 0, that is we re- 
produce all the desired constraints ([TP]). ([TTj). (flUl) . (f2T))) . 
The constraints Ti, T 2 have vanishing Poisson brackets 
with all the constraints. The remaining constraints obey 
the Poisson-bracket algebra 

{T 3 , T 4 } = -4T 5 , {T 3 , T 5 } = -2T 3 + 2a 3 , 

{T 4 ,T 5 } = 2T 4 -2a 4 . (23) 



If we take the combination 

T 3 = T 3 + -^T 4 , 
a 4 

the algebra acquires the form 



(24) 



{T 3 , T 4 } = -4T 5 , {T 3 , T 5 } = -2T 3 + 2^T 4 , 

a 4 

{T 4 ,T 5 } = 2T 4 - 2a 4 . (25) 

The only bracket which does not vanish on the constraint 
surface is {T 4 , T5}. According the Dirac terminology [13- 
15], we have the first-class constraints ([M]), ([H]), (|2T))) . 
and the second-class pair (fTTj) . The presence of the first- 
class constraints indicates that we are dealing with a 
theory invariant under a three-parameter group of local 
(gauge) symmetries, which will be found below. 

The constraints ([24]). (ITU can be taken into account 
by transition from the Poisson to the Dirac bracket, and 
after that they can be omitted from consideration [13- 
15]. Since the Dirac brackets are constructed with use 
of 5*0(2, 3) -invariants, the Dirac brackets of the quan- 
tities J coincide with the Poisson one ([8]). Hence we 
quantize the model according Eq. (fT5|) . The operator 
produced by the first-class constraint C|19[) is imposed on 
the state vector. This gives the Dirac equation. In the 
result, canonical quantization of the model leads to the 
desired quantum picture. 

The Lagrangian formulation can be restored from 
Eq. (j21[) . To achieve this, we note that the conju- 
gate momenta enter into the Hamiltonian in the form 
\P a G a ^Pfi, where P a = tt v , 7r 5 ) and the "metric" 
G a ^(w A1 ,w 5 , ei, e 2 , e 3 ) is given by 9 x 9 -matrix 



e 2 u} 5 rj 



I -e 2 c \ 



-e 2 w u 







1 




e 3 ?7 

1 


1 •■• 

1 











-ea / 



Here 77 is the Minkowski metric. To find the Lagrangian, 
we solve the Hamiltonian equations for the position vari- 
ables Q a = (x^w^w 5 ), Q a = G af3 P p , with respect to 
P a . It gives P a = G a pQP, where G a p is the inverse met- 
ric. We substitute these P a back into Eq. (|2ip . obtaining 



the Lagrangian 
1 



L = ^G a0 Q a Q? 



-e 4 w lo a - -e b a b . 



(26) 



It has been denoted (a±, a 2 , a 3 , a 4 ) = (m 2 c 2 ,mch, a 3 , o 4 ). 
The kinetic term looks like those of a free particle moving 
on the curved nine-dimensional space with the metric 
G a p. 

We introduce the abbreviation 



e 2 



(27) 



where = lu 5 lu^ — lo^lo 5 is the configuration-space 
counterpart of the angular momentum J 5ai . Then mani- 
fest form of the Lagrangian (j2GI) is 



B 



1 

f J 3 



{Dx»f - ^(Dx^) 2 



A ■ A 



(28) 



It has been denoted B = eie 3 -e 2 (oj 5 ) 2 , A = B+e\{u^) 2 . 
We have verified that hamiltonization of the Lagrangian 
leads back to the Hamiltonian action (l2~Tj) . 



IV. CLASSICAL EQUATIONS OF MOTION 
AND THE LACK OF ZITTERBEWEGUNG 



In this section we show that the status of Zitterbewe- 
gung phenomenon in relativistic quantum mechanics is 
not as clear as is commonly believed. 

Besides the constraints discussed above, the action 
(j2"Tj) implies the following equations (we use the notation 
(pu>) = p^ujf, ) 



&a — ^ea? ^ea — 0, G. — 1, 2, 3] 

C4 — — e 3 , 7T e4 = U, A e4 = — A e3 ; 
a 4 a 4 



(29) 



-e 2 J 5 ^, p"=0; 



tu^ = e 3 vr M + e 2 wV\ tt' 1 = e 2 7rV - — 



lj 5 = e 3 7r 5 + e 2 (puj), 7r 5 = e 2 (pn) -e 3 u; 5 . (30) 

a 4 



As a consequence, we obtain (assuming the gauge condi- 
tions e a = const) 'x^ = —h^J^ v p u , that is x^ -variable 
is under Zitterbewegung in non-interacting theory. 

Since we deal with a constrained theory, our first task 
is to specify the physical-sector variables. We note that 
Eqs. (j29|) . (j30j) do not determine the Lagrangian multi- 
pliers A a (r), a = 1,2,3. Eqs. (j2"9")l then imply that e a , 
a = 1,2,3, remain an arbitrary functions as well. The 
arbitrary functions enter into solutions to equations of 
motion for the variables u> A , ir A . Hence, except 
all the variables has ambiguous dynamics. According to 
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the general theory [13-15], variables with ambiguous dy- 
namics do not represent the observable quantities. In 
particular, the coordinate x^, which corresponds to the 
center-of-charge position operator of the Dirac equation, 
and experiences Zitterbewegung, is an unobservable quan- 
tity in our model. 

Let us compute the total number of physical degrees of 
freedom. Omitting the auxiliary variables and the corre- 
sponding constraints, we have 18 phase-space variables 
x M , oj a , tta subject to the constraints (fTO]) . ((TTJ) , 
([T9| . ([20]). Taking into account that each second-class 
constraint rules out one variable, whereas each first-class 
constraint rules out two variables, the number of physi- 
cal degrees of freedom is 18 — (2 + 2x3) = 10. To con- 
struct the unambiguous position variable, we note that 
the quantity 

a" = x» + ^J^Pu, (31) 
2p z 

obeys x 1 = ep^, where e = e\ + Besides, we 

know = 0. Since these equations resemble those for 
a spinless relativistic particle, the remaining ambiguity 
due to e has the well-known interpretation, being re- 
lated with reparametrization invariance of the theory. In 
accordance with this, we assume that x^ir) represent 
the parametric equations of the trajectory x l (t). Us- 
ing the identity dA ^ = c ±o^ > we conclude that the 
reparametrization- invariant variable x l (i) has determin- 
istic evolution: = cfn = . In the absence of inter- 
action, it moves along a straight line. We also notice that 
represents the center-of-mass (Pryce-Newton-Wigner) 
coordinate [3, 4], and p^ represents its mechanical mo- 
mentum. Hence the mass-shell condition (|2Tfl) guarantees 
that the i l -particle cannot exceed the speed of light. 

As the classical four-dimensional spin vector, we take 
the Pauli-Lubanski vector 5 M = -z^^PvJap- It has no 
precession in the free theory, — 0. In the rest frame 
p^ = (mc, 0, 0, 0), it reduces to the three-dimensional ro- 
tation generator, 5° = 0, S 1 = \mce^ k Sjk, as is ex- 
pected in the non-relativistic limit. Hence the ten vari- 
ables x l (t), p % {t), S^it) can be taken as the physical vari- 
ables. 

We have specified the physical sector from analysis of 
equations of motion. The more traditional way to do 
this consists of analysis of local symmetries of the for- 
mulation. The action (|21[) is invariant under a three- 
parameter group of local symmetries. One of them is the 
reparametrization symmetry which we take in the form 

S a Z = a{Z, H}, S a e a = (ae a ). (32) 

Here Z — (x^,p^,uj a ,tt a ). This form of reparametriza- 
tion symmetry can be obtained as follows. It is sufficient 
to consider the spinless particle action (1221) . The stan- 
dard form of reparametrization symmetry reads Sx^ = 
ax 1 *, Sp^ = ap^ , Se — (ae) (we have omitted the trans- 
formations 5X e — (Se)' and <57r e = which are not 



relevant for our discussion). We rewrite the symme- 
try in equivalent form, without derivatives acting on x 
and p. Every Hamiltonian action has trivial symme- 
tries [20], in this case they are 8x^ L — e(x^ — {x^,H}), 
Sp^ = e(p^ — {p^,H}), where H stands for the com- 
plete Hamiltonian. Taking the combination 5 a + S e with 
e = —a, the reparametrization invariance acquires the 
form Sx^ — a{x^,H}, Sp^ = alp 11 ,!!}), Se = (ae). 

For our case (f2~Tj) . the manifest form of the 
reparametrization symmetry is 

fe" = a( ei p" + \e 2 J 5 »), Sp» = 0, 5e a = (ae a ), 
Slo" = a(e 3 7r' 1 + e 2 u°p^), Sw" = a(e 2 7rV - — e 3 ^), 
Slj 5 = a(e37r 5 + ezOpu))), Sir 5 — a(e 2 (pTr) -e^uj 5 ). 

Besides, there are two more symmetries with the local 
parameters /3(r), 7(7-) 

8pxP = /V\ Spex = $; (33) 

S lU j A = jn A , (5 7 e 3 = 7. (34) 

All the initial variables, except p M , are not gauge in- 
variant. The center of mass x^ turns out to be in- 
variant with respect to two symmetries: <5 7 i M = 0, 

SpH* = 0, where ^ = Sp + S a , a = - ^^^ji) ■ 
As expected, i M is affected only by reparametrizations: 
S a x^ = a(ei + C2 2p2 )P^- The spin vector S 1 * is invari- 
ant under all the transformations. So we have confirmed 
our previous result: x M , p^ and can be taken as the 
physical-sector variables. 

We finish with a preliminary comment on interac- 
tion with an external electromagnetic field. The clas- 
sical constraints that produce Eqs. ((T|), ((2]) are T 2 = 
(p M + fA M )J 5 ^ + mch = 0, Ti = (pP + f^) 2 + 
■ 2 \F^vJ llv + m 2 c 2 — 0. Their Poisson bracket reads 
{T 2 ,Ti} = -T^d a F^J 5a J^. For the homogeneous elec- 
tric and magnetic fields, the constraints form the first- 
class system, {T 2 ,Ti} = 0. Hence the interaction does 
not break the local symmetries presented in our model. 
In the general case, the breaking is in a sense "soft", i.e., 
proportional to ft 2 . Hence, to construct an interaction 
with an arbitrary field, one can start the iteration proce- 
dure, adding the non-minimal-interaction terms of order 
h 2 or more to the constraints T%, T 2 . 

V. CONCLUSION 

In this work we have constructed a semiclassical model 
(f2"Tj) . (|2"5)) for description of the relativistic spin and 
showed its consistency both on the classical and on the 
quantum level. Canonical quantization of the model 
leads to the Dirac equation. As we could expect for the 
relativistic spinning particle, the physical sector of the 
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model is composed by the position variable x^ (l3l"j) and 
the spin vector S^. In the absence of interaction, they 
obey the free equations x" — 0, = 0. 

We have presented a simple semiclassical argument 
that prohibits the relativistic Zitterbewegung. Roughly 
speaking, the argument is as follows. The Dirac equa- 
tion (pj implies the Klein- Gordon one ([2]). In contrast, 
in the classical theory the corresponding constraint (fT9|) 
does not imply the mass-shell constraint p 2 + m 2 c 2 = 0. 
To obtain a consistent picture, we are forced to construct 
the semiclassical model that produces both constraints. 
In turn, the presence of independent constraints implies 
that we deal with a theory with the local symmetries 
(f32 j) - (j34|) . Physical quantities are those invariant under 



the local symmetries. Our observation is that the classi- 
cal variable x^, that corresponds to the center-of-charge 
position operator of the Dirac theory, and experiences 
Zitterbewegung, is not invariant. Similarly to the poten- 
tial of electromagnetic field, x^ is an unobservable 
quantity. 



VI. ACKNOWLEDGMENTS 

This work has been supported by the Brazilian foun- 
dation FAPEMIG. 



[1] E. Schrodinger, Sitzunger. Preuss. Akad. Wiss. Phys.- 

Math. Kl. 24 (1930) 418. 
[2] P. A. M. Dirac, The Principles of Quantum Mechanics 

(Clarendon Press, Oxford, 1958) p. 257. 
[3] M. H. L. Pryce, Proc. Roy. Soc. A 195 (1948) 62. 
[4] T. D. Newton and E. P. Wigner, Rev. Mod. Phys. 21 

(1949) 400. 

[5] L. L. Foldy and S. A. Wouthuysen, Phys. Rev. 78 (1950) 
29. 

[6] J. Frenkel, Z. fur Physik 37, (1926) 243. 

[7] V. Bargmann, L. Michel and V.L. Telegdi, Phys. Rev. 

Lett. 2 (1959) 435. 
[8] A. J. Hanson and T. Regge, Ann. Phys. 87 (1974) 498. 
[9] F. A. Berezin and M. S. Marinov, JETP Lett 21 (1975) 

320; Ann. Phys. 104 (1977) 336. 
[10] A. O. Barut and N. Zanghi, Phys. Rev. Lett. 52 (1984) 

2009. 

[11] H. C. Corben. Classical and Quantum Theories of Spin- 
ning Particles (Holden-Day, San Francisco, 1968); I. 
M. Ternov, Introduction to Physics of Spin of Elemen- 
tary Particles (Moscow State University Publ. House, 
Moscow, 1997) (in Russian). 

[12] B. P. Kosyakov and V. V. Nesterenko, Phys. Lett. B 
384 (1996) 70; A. Berard and H. Mohrbach, Phys. Lett. 
A 352 (2006) 190; N. Kudryashova and Y. N. Obukhov, 
Phys. Lett. A 374 (2010) 3801. 



[13] P. A. M. Dirac, Lectures on Quantum Mechanics (Yeshiva 
University, New York, 1964). 

[14] D. M. Gitman and I. V. Tyutin, Quantization of Fields 
with Constraints (Springer- Verlag, Berlin, 1990). 

[15] A. A. Deriglazov, Classical Mechanics, Hamiltonian and 
Lagrangian Formalism (Springer- Verlag, Berlin Heidel- 
berg, 2010). 

[16] A. A. Deriglazov, Phys. Lett. B 626 (2005) 243. 

[17] A. A. Deriglazov, Mod. Phys. Lett. A 25 (2010) 2769. 

[18] J. J. Sakurai, Advanced Quantum Mechanics, (Addison- 
Wesley Publishing Company Inc., 1967) p. 139. 

[19] G. Salesi and E. Recami, Phys. Lett. A 267 (2000) 219; 
W. Zawadzki and T. M. Rusin, Phys. Lett. A 374 (2010) 
3533; P. Kurzynski, Phys. Lett. A 372 (2008) 6125; O. 
Roslyak, G. Gumbsa and D. Huang, Phys. Lett. A 374 
(2010) 4061; S. Sasabe and K. Tsuchiya, Phys. Lett. A 
372 (2008) 381; 

[20] A. A. Deriglazov and K. E. Evdokimov, Int. J. Mod. 
Phys. A 15 (2000) 4045; A. A. Deriglazov, J. Phys. A 
40 (2007) 11083; J. Math. Phys. 50 (2009) 012907-1-15. 

[21] A. A. Deriglazov, Ann. Phys. 327 (2012) 398; Arxiv: 
1107.0273. 

[22] F. Zahringer, G. Kirchmair, R. Gerritsma, E. Solano, R. 
Blatt, and C. F. Roos, Nature (London) 463 (2010) 68. 



